Structural changes in quasi- ID many-electron systems: from linear to zig-zag and 

beyond 
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Many-electron systems confined to a quasi- ID geometry by a cylindrical distribution of positive 
charge have been investigated by density functional computations in the unrestricted local spin 
density approximation. Our investigations have been focused on the low density regime, in which 
electrons are localised. The results reveal a wide variety of different charge and spin configurations, 
including linear and zig-zag chains, single and double-strand helices, and twisted chains of dimers. 
The spin-spin coupling turns from weakly anti-ferromagnetic at relatively high density, to weakly 
ferromagnetic at the lowest densities considered in our computations. The stability of linear chains 
of localised charge has been investigated by analysing the radial dependence of the self-consistent 
potential and by computing the dispersion relation of low-energy harmonic excitations. 



I. INTRODUCTION 

Systems made of ionised atoms confined into Pen- 
ning ll| or Paul [2!] traps provide interesting prototypes 
of low-dimensional many-particle systems. Confinement 
is enforced by the application of static (Penning) and/or 
time-dependent electromagnetic fields (Paul) and tem- 
perature can be reduced to the mK range by laser cool- 
ing [31 . Systems composed of atoms ranging in number 
from a few tens to several thousand are routinely made 
and their properties analysed by a broad range of spec- 
troscopic techniques [J]. 

Fine control of the ion density pi and temperature T 
allows to experimentally probe a wide range of inter- 
particle couplings, whose strength is measured by the 
ratio r between the potential and kinetic energy of the 
ions. At sufficiently high coupling (F ~ 180, see Ref. 0), 
ions condense into a regular lattice, fulfilling early pre- 
dictions generally attributed to Wigner 0] . Such a tran- 
sition has been observed in trapped clouds made of sin- 
gle [7| or binary J8| ion species, providing an intriguing 
view of ordered Coulomb systems at low temperature and 
high couplings [3,[i3|- Upon changing pi, T or the shape 
and strength of the applied fields, these systems undergo 
a sequence of characteristic structural changes such as 
order-disorder ill' and isomerisation transitions [ij] . 

Recently, experiments |13l , ll4| and computational stud- 
ies [l3 have focused on one specific phase change taking 
place in nearly ID trapped-ion systems, transforming lin- 
ear chains into zig-zag configurations [1(| . Simple consid- 
eration of the forces active in and on the system suggests 
that such a transition arises from the competition be- 
tween the interaction with the external field (increasing 
upon the transition) and the electron-electron repulsion, 
which decreases because of the larger nearest neighbour 
distance in zig-zag chains. The order of the transition 
is not precisely known and, strictly speaking, not even 
well defined for the finite samples probed in experiments. 
However, computations for extended systems with peri- 
odic boundary conditions suggest that the transition is 
continuous and second order, with a discontinuity in the 



second derivative of the ground state energy with respect 
to the ID ion density |15| . 

The linear to zig-zag transition has a number of dif- 
ferent implications on the system properties. In the 
case of ions carrying a magnetic (spin) moment, for in- 
stance, changing the amplitude of the zig-zag modulation 
changes the number and relative distance of the ions' 
neighbours, and, by varying the relative weight of first 
and second nearest neighbour interactions, it provides a 
way to tune the spin-spin coupling [13] ■ At the same 
time, the zig-zag transition doubles the unit cell of the 
system, thus changing even qualitatively the vibrational 
spectrum of the chain. Interestingly, the dimerisation 
that often accompanies the zig-zag transition could open 
the way to the Bose-Einstein condensation of ions whose 
spin is half integer. Finally, the link between the struc- 
tural transitions seen in experiments and the behaviour of 
quantum correlations among the trapped ions has been 
established by a theoretical study based on techniques 
typical of continuous- variable systems [lal ■ 

Besides providing an appealing playground to inves- 
tigate the interplay between dimensionality and many- 
particle effects, quasi-lD atomic plasmas are actively in- 
vestigated in view of applications in metrology [19| and 
in quantum information technology [20]. Moreover, low 
temperature ions confined in a quasi-lD ion trap have 
been proposed as models for the analogic simulation of 
many-particle systems |2l| - |23| . In most of these applica- 
tions, quantum mechanics plays an important role. This 
observation has motivated us to study quasi- ID systems 
made of electrons, whose light mass amplifies the quan- 
tum mechanical effects. To the best of our knowledge, 
low-dimensional many-electron plasmas confined at low 
temperature into electromagnetic traps have not been 
made and characterised in experiments. However, it 
should be mentioned that there is a considerable interest 
in scaling up the already experimentally demonstrated 
ability to trap and coherently control a single electron in 
a Penning trap to a genuine quantum many-body config- 
uration. Roadmaps towards the achievement of such a 
situation have been detailed in Ref. [2g]. It is worth re- 
minding that closely related systems, consisting of mobile 



electrons in conducting nano-wires can be prepared by a 
variety of methods, including the controlled doping of 
semiconductingnanostructures [2J, [23| and conducting 
polymers |27h29| . Moreover, new fabrication methods 
are being developed based on the doping or the electro- 
static biasing of carbon nanotubes [30| • 

The properties of electrons in these systems are often 
described using the one-component plasma picture [3l|, 
thus representing electrons as independent particles mov- 
ing in a fixed external potential. We resort to the simplest 
version of this approach, based on the jellium model [S^] , 
in which the external potential confining the electrons 
is due to their Coulomb interaction with a cylindrical 
background of positive charge, whose density is constant 
(pb) within a pre-defined volume and zero outside. In 
order to approach the conditions of interest for charged 
particles in a trap, we consider the limit of very thin 
wires with a large aspect ratio between the length L^ 
and the radius Ri, of the background charge distribution 
(27.7 < Lz/Rb < 39.2). In our computation, many-body 
effects are accounted for by resortin g to the simple local- 
spin density (LSD) approximation |33{ to density func- 
tional (DF) theory 34]. Single electron orbitals are ex- 
panded on a large basis of plane waves, and the ground 
state energy and density are determined by direct min- 
imisation, without any symmetry restrictions. We focus 
our attention on the low density, high correlation regime, 
where the Wigncr crystal is the stable phase, and we 
consider various combinations of spin-up and spin-down 
populations. 

Our calculations demonstrate the existence of linear 
and zig-zag chains, stabilised by different combinations of 
the ID electron density, spin configuration and shape of 
the confining potential. Furthermore, our results display 
other unforseen structures, never considered or found so 
far for classical ion systems. Dimers, already suggested 
long ago for quantum spin chains j35| , appear at low den- 
sity. At intermediate densities, we find new geometries, 
such as helices, and even double helices. Hints of these 
exotic geometries were already given by calculations for 
larger wires at much lower density |36j . 

In addition to this basic information on the ground 
state density and spin configuration, our results provide 
a wealth of new data on the density of states, the elec- 
tric conductivity and the vibrational modes of nearly ID 
electron systems. In particular, we give strong evidence 
of a second-order nature of the linear-to-zig-zag transi- 
tion, a result that appears to be in line with the findings 
in Ref. [15]. 

The paper is organised as follows. The model and 
the computational method are defined and briefly dis- 
cussed in Section ini The computational results for the 
ground state density and spin distribution are described 
in Sec. IIII Al while the computation of phonon-like exci- 
tations is reported in Sec. IIIIBl A summary and a brief 
outline of promising new directions are given in Sec. IIVI 



II. MODEL AND METHODS 



Computations have been carried out for systems of 
N = Nup + Ndn electrons, neutralised by a cylindri- 
cal background of positive charge, whose axis is paral- 
lel to the z direction. In what follows, the background 
density pb is expressed in terms of the Wigner-Seitz ra- 
dius Ts through the relation pb — 2>/^Tir^g. Here Nup 
and Ndn are the number of spin-up and spin-down elec- 
trons, respectively. The length L^ and radius Rb of the 
cylindrical background satisfy the neutrality condition 
irR^LzPb = N. Moreover, the number of electrons per 
unit length of the wire is pim = nR^pb. Atomic units are 
used throughout the paper, and cylindrical coordinates 
(r, 0, z) are implicitly assumed in our equations and de- 
scription of the results. 

The basic cylindrical segment described above is peri- 
odically replicated in the direction parallel to the z axis 
with periodicity L^, thus representing an extended wire 
along such direction. Due to our choice of plane waves as 
basis functions (see below) , we periodically replicated our 
sample also in the xy plane. For the sake of simplicity, 
we adopt the same periodicity L^ in all three directions. 

The ground state energy and density are computed 
within the Kohn-Sham (KS) formulation of density func- 
tional theory (DFT), in which electrons occupy single- 
particle KS states {V'i;* = l,...,N}. The density p(r) 
and spin polarisation m(r) are given by 
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The ground state energy and density are determined by 
minimising the KS energy functional 
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where Vext{r) is the Coulomb potential of the positive 
charge distribution. Here Uxc[p\ is given by the local 
spin density approximation: 



Uxc[p]^ /p(r)exc(p(r);C(r))dr, 



(4) 



where exc(p(r); C(r)) is the exchange-correlation energy 
per electron [3j] of the homogeneous electron gas at 
the local density p(r) and local spin polarisation C(r) = 
m{v)/p{v). 

Kohn-Sham orbitals are expanded on a basis of plane 
waves whose periodicity matches the cubic periodicity of 



the simulation cell: 
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where each G is a reciprocal lattice vector of the cubic 
simulation cell. The basis includes all plane waves whose 
G vector satisfies G^ < Ecut with Ecut a suitable kinetic 
energy cut-off. A plane wave basis set of cut-off 2Ekin is 
used to represent the electron density and the external 
potential Vexti^)- The Fourier expansion of the latter is 
easily obtained by using Poisson's equation, and consid- 
ering that the Fourier transform of the positive charge 
density is given by 
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if G = 0, 
otherwise, 



(6) 



where V = L^Lz, and Ji is the cylindrical Bessel function 
of order one. As implicit in our notation in Eq. ([Ij-Q, 
the Brillouin zone defined by the periodicity of the simu- 
lated system is sampled at the T point only. The size of 
the systems we simulate is such that this approximation 
does not introduce any sizable error. We also verified 
that the density overlap and the spurious interactions 
across the xy plane are negligible. 

The optimisation of the Eks[p] functional is carried 
out by direct minimisation, i.e. by considering Eks[p\ 

as an algebraic function of the Fourier coefficients Cq , 
and using standard minimisation routines |37| . In do- 
ing so, we use the following expression for the derivative 
of the Kohn-Sham energy functional with respect to the 
expansion coefficients {cq } 
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The last equality in Eq. ([7]) implicitly defines the Kohn- 
Sham Hamiltonian HkS: which can also be written as 



Hks = -Jv^ 



VKsir), 



(8) 



where VR-s(r) is the self-consistent KS potential. 



III. RESULTS OF THE SIMULATIONS 



A. The ground state density distribution 

Computations have been carried out for systems of 
16 to 32 electrons, neutralised by a cylindrical back- 
ground of aspect ratio 27.7 < Lz/Rb < 39.2. Differ- 
ent total spin polarisations have been considered, from 
Nup = Ndn = N/2 to Nup = N, Ndn = 0. Our compu- 
tations, however, are spin unrestricted, and even in the 
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FIG. 1: (Color online) Density iso-surface p = 2pi, for samples 
oiN = 32 electrons, (a) r, = 30; (b) r^ = 32; (c) r^ = 35; (d) 
Vs = 40. The vertical arrow marks the position of a defective, 
unpaired charge blob (see text). 



Nup = Ndn case, spin polarisation can arise locally driven 
by exchange interactions. 

We consider systems of fairly low density, correspond- 
ing to 25 < Ts < 40. Previous studies on the homoge- 
neous 3D electron gas, carried out using the same ap- 
proach as ours, have shown that DFT-LSD predicts the 
onset of spin polarisation and charge localisation to take 
place at r^ ^ 25 [33, Is^ , while electrons appear to be 
well localised at rj > 30. 

Our approach provides directly the ground state en- 
ergy and density distribution for any given spin popu- 
lation as a function of size and shape of the background 
density. The results for the total energy, however, display 
fairly predictable trends as a function of the background 
density and linear density pun- For this reason, in full 
analogy with what is done in the classical case of atomic 
ions, we focus our discussion on the dependence of the 
density and spin distribution on the model parameters. 

The linear to zig-zag transition in experimental quasi- 
ID ion systems is usually triggered by changing the geo- 
metric shape of the external potential which confines the 
ions. As expected, stiff potentials that tightly confine 
ions in the radial direction, favour linear configurations, 
while soft potentials give rise to zig-zag chains. In our 
model, the curvature of the external potential along the 
radial direction is directly related, via Poisson's equation, 
to the background density, and decreases with decreasing 
Pb- Therefore, we first present the results of computa- 
tions for samples of increasing r^ (decreasing pb), keep- 
ing fixed the number of electrons, the total spin, and the 
aspect ratio L^/Rb- We consider, at first, fully spin po- 
larised systems, which arguably represent the simplest 
case. 

The electron distribution is characterised by plotting 
density iso-surfaces, which in the homogeneous (3D) elec- 
tron gas case show a regular pattern of charge blobs for 
rs > 30 (see Ref. ^). The results for Nup = N ^ 32, 
Ndn = at r^ = 30, 32, 35, and 40 are shown in Fig.[T](a) 



- (d). For all these systems the periodicity of the simu- 
lation cell along the wire axis is L^ — 32rs, and the ratio 
L^l Rb is 27.7. As already stated, the same periodicity is 
used along the other two directions, resulting in a large 
super-cell {V = 32768rf), and a fairly high number of 
plane waves {upw '^ 130000) in the expansion of the KS 
orbitals. 

The Ts — 30 ground state density displays a clear he- 
lix geometry, with 17 full periods within the simulation 
cell. The helix period is close but not equal to the length 
A — l/{2kp) expected on the basis of known singularities 
of the response function [401 , suggesting that at such low 
densities localisation cannot be quantitatively described 
in terms of linear response. In the equation above, kp 
is the Fermi wave vector in the z direction, evaluated by 
computing KS bands in the ID Brillouin zone, and as- 
suming a cylindrically symmetric and translationally in- 
variant charge density along the wire 31]. By comparing 
the number of particles and the number of helical turns in 
the simulation cell we see that slightly less than two elec- 
trons are accommodated in each turn. The non-integer 
number of electrons per turn suggests that the system 
might not be a closed-shell configuration. Therefore, the 
addition or the subtraction of electrons, while adjusting 
the background parameters to keep the system neutral, 
could enhance the ground state stability. We remark 
again that at Ts ~ 30, DFT-LSD for the homogeneous 
electron gas gives a well localised ground state charge 
distribution at all spin polarisations. The helix config- 
uration, therefore, results exclusively from the quasi-lD 
confining potential. 

The discretisation of the charge density into localised 
blobs becomes apparent again already at r^ — 32. Lo- 
calisation, however, is still incomplete, and the charge 
distribution gives rise to a crankshaft-like structure par- 
allel to the wire axis. The crankshaft harms are marked 
by pairs of elongated and partially overlapping blobs of 
electronic charge. According to the results of the min- 
imisation process, the number of blobs in the simulation 
cell is 47, i.e., significantly higher than the number of 
electrons in the system. 

The partition of charge into blobs becomes progres- 
sively more marked with decreasing density, and is com- 
plete at Ts = 35 (see Fig. [T] (c)). At r^ = 40 (Fig. [1] 
(d)) the charge distribution can be described as a line 
of charge dimers whose direction displays a rather com- 
plex pattern in space. A closer analysis, however, reveals 
that this pattern can be described as due to a full rota- 
tion of the dimer bond about a direction perpendicular 
to the wire axis. The dimer rotation is refiected in the z 
dependence of the in-plane averaged density, defined as 
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and shown in Fig.[2]for Vg = 40. Each of the density peaks 
in the rapidly oscillating part of this plot (0 < z < lOrs) 
marks the position of dimers perpendicular (or nearly 
perpendicular) to the wire axis. The nearly constant 
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FIG. 2: (Color online) Upper panel: planar average of the 
electron density p{z) (see the definition in Sec. illlll . Lower 
panel: planar average of the kinetic energy density ek[z) di- 
vided by the corresponding density p{z). ek{z) is defined in 
analogy with the planar density p{z). 



portion at z '^ br^ corresponds to a few dimers nearly 
parallel to the z axis. A similar behaviour is displayed 
by the in-plane average of the electron kinetic energy, 
shown in the upper panel of Fig. [2] for a comparison. 

The number of blobs remains constant at 47 with de- 
creasing density from Ts = 32 to r^ = 40, thus preventing 
the identification of blobs with single electrons. Clearly, 
the odd number of blobs implies that not all of them form 
dimers. One defective blob, in fact, is easily located and 
is identified by the vertical arrow in Fig. [1] (d). Apart 
from such an isolated defect, the structure of the fully 
polarised system at r^ = 40 is fairly regular, with an 
intra-dimer separation almost exactly equal to Tj and a 
dimer-dimer separation of I.4rs. 

Our discussion has been focused, so far, on the de- 
pendence of the electronic structure on the curvature of 
the external potential. The shape of orbitals, however, is 
more directly related to the spatial extension of the self- 
consistent KS potential, which is very sensitive to the ID 
density pun- To explore this dependence, we carried out 
computations for systems of the same length Lz , but dif- 
ferent background radius i?(,, corresponding to systems 
having a number of electrons between 16 to 32. The 
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FIG. 3; (Color online) Density iso-surface p — 2ph for fully 
spin polarised samples of equal length L^, and different Rb 
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data for the density distribution at Ts = 30 and full spin 
polarisation {Nup — N) are displayed in Fig. [3l which 
shows that with increasing pun the distribution of charge 
blobs goes from linear towards more complex shapes, the 
transition taking place in between TV = Nup = 24 (lin- 
ear) and iV = Nup = 26 (bent). More in detail, the 
TV = 16 and A^ = 18 samples show a cylindrical charge 
distribution, iV = 20 to 24 correspond to linear chains of 
blobs, elongated at first, and then progressively rounded 
with increasing N. At N = 26 the structure is a zig-zag 
chain, transforming into an helix above A^ = 28. The 
sequence of ground state structures as a function of N is 
similar for Ts = 40, with the exception that, at high N^ 
the helix apparent in the rg = 30 data is replaced by a 
string of dimers at r^ =40. The transition from linear 
to non-linear configurations, in particular, takes place at 
the same size {N — 24) at r^ = 30 and r^ — 40. 

As already suggested, the origin of these changes can 
be traced back to the Kohn-Sham potential, whose de- 
pendence on the radial distance is shown in Fig. 2] The 
progressive widening and softening of VKsi'>') with in- 
creasing Tg is apparent in this figure and can be quantified 
by fitting Vfs-5'(r) at short r with the sum of a quadratic 
and a quartic term. The fitting coefficients, given in the 
inset of Fig. HI show that the quadratic term becomes 
rapidly less important than the quartic one with increas- 
ing A^, until it vanishes for A^ ~ 24. At the same time, 
the coefficient of the quartic term increases with increas- 
ing A^, thus preserving the overall stability of the system, 
at the expense of the linearity of the chain. 

Comparison of the charge distribution with the in- 
formation given by VR-5(r) suggests that the cylindrical 
charge distribution for A^ = 16 — 18, as well as the elon- 
gated shape of blobs seen at A^ = 20 is apparently due 
to the squeezing effects of a narrow harmonic potential. 
For A^ > 24, when the configuration is bent, confinement 
is exclusively due to the quartic term, while the negative 
quadratic term gives rise to an off-centre minimum. At 



FIG. 4: (Color online) Dependence of the Kohn-Sham poten- 
tial on the radial distance r for fully spin polarised samples 
of equal length L^ and different Rt at r^ — 30. (and thus 
of different values of the ID density pun)- In order of in- 
creasing Vks{0), the curves correspond to: Ai' = 16, A'' = 20, 
iV = 24, iV = 28, A'' = 32. Inset: iV-dependence of the 
quadratic (02) and quartic (04) coefficients in the polynomial 
fit of VKs{r) along the radial direction. The a2 coefficients 
have been rescaled (multiplied by 1000) to plot a2 and 04 on 
the same scale. 
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FIG. 5: (Color online) Spin polarisation iso-surface m{r) = 
2 pb at background density corresponding to Vs = 30 for sam- 
ples of iV = 32 electrons and different spin populations, (a) 
Nup = 32, A^dn = 0; (b) N^p = 31, Nan = 1; (c) N^p = 16, 
Ndn ~ 16. Blue (dark) surface: spin up electrons. Yellow 
(light) surface: spin-down electrons. 



the size of the transition (A^ — 24), the short range por- 
tion of VR-s(r) is very flat, and is well represented by a 
r^ term. Further insight into the stability of the ground 
state structure found by our minimisations is obtained 
by changing the net spin of the sample. 

Energy differences among configurations of the same 
size and Vs but different Nup, Ndn tend to be small at 
the densities considered in our study. Nevertheless the re- 



suits of our DFT-LSD computations shown that the spin- 
spin couphng turns from anti-ferromagnetic at r^ = 25 
to ferromagnetic for r^ > 32. The evolution of the mag- 
netic structure upon changing the ratio of N^p and Ndn 
at fixed background density is illustrated in Fig. [S] dis- 
playing magnetisation iso-surfaces for samples of iV = 32 
electrons at r^ = 30. Starting from the helical struc- 
ture of the fully polarised case (Fig. [5] (a)), reversing one 
spin in the r^ = 30 case gives rise to a localised charge 
and spin blob that breaks the continuity of the helix (see 
Fig. [5] (b)). Even more striking and, at the same time, 
more significant, is the result obtained by reversing half 
of the 32 electron spins. In the r^ = 30 case, this breaks 
the single helix into a double helix, whose two strands 
have opposite spin (see Fig.[5](c)). 

Reversing one spin on the Vs =40, fully polarised sam- 
ple localises the reversed spin on the single blob not fit- 
ting into the dimer pattern described above (see Fig. [6]) . 
This observation confirms our identification of this un- 
paired blob as a defect, whose stability is intrinsically 
lower than that of the dimerised blobs. Reversing now 
half of the spins in the Vg = 40, TV = 32 sample re- 
sults into a configuration somewhat similar to that of the 
Ts = 30, Nup = Ndn case. At this low density, however, a 
sizable amount of disorder makes the identification of the 
underlying double helix pattern more difficult. The en- 
hancement of disorder is probably due to the further de- 
crease of the exchange coupling with decreasing density. 
A comprehensive view of the dependence of density and 
magnetic structures on pb (or equivalently rg ) in globally 
spin-compensated samples is given in Fig. [T) 

The density of states computed from the Kohn-Sham 
eigenvalues in all cases consists of a few clearly identi- 
fiable bands, some of them overlapping. Especially at 
low density, some of the bands are disjoint. The results 
for fully spin-polarised systems shown in Fig. [8] confirms 
our anticipation that these samples are open-shell sys- 
tems, whose Fermi energy falls in the middle of a band. 
The quantitative picture emerging from computations, 
however, depends on size and spin, and in the case of 
the spin-compensated samples considered in our study 
{Nup — Ndn — 16) the Fermi energy falls into a well 
defined and fairly wide gap for rs > 30. 

Conductivity has been computed using the Kubo- 
Greenwood formula, upon computing a fairly large num- 
ber of emp ty states using the method briefiy discussed 
in Ref. (Sg]. The results (not shown) refiect the features 
seen in the DOS. The low frequency conductivity is high- 
est for ferromagnetic samples, apparently because of their 
open-shell character. Moreover, in most cases we find 



FIG. 6: (Color online) Spin polarisation iso-surface m{r) = 
2p6 at background density corresponding to rs — 40. System 
of 32 electrons, Nup = 31, Ndn = 1- Blue (dark) surface: spin 
up electrons. Yellow (light) surface: spin-down electron. 



non-negligible conductivity down to fairly low frequency, 
due to the presence of a few structural defects in the 
ground state electron distribution. 



B. Low^-energy phonon-like excitations 

The subdivision of charge into weakly overlapping 
blobs motivates us to investigate the possibility of ob- 
serving phonon-like excitations, corresponding to small 
displacements of the centre of mass of individual blobs. 
This analysis, however, appears to be meaningful only for 
the cases where the number of blobs corresponds exactly 
to the number of electrons in the system, in such a way 
that we can identify blobs with single electrons. Our dis- 
cussion below concerns one of such cases, corresponding 
to Nup = A^ = 24, r^ = 40. This system is in fact of par- 
ticular interest, since it marks the transition from linear 
to zig-zag configurations, and phonons are expected to 
reflect the impending change in the charge distribution. 

To estimate phonon frequencies, we approximate 
Kohn-Sham orbitals with single Gaussians, centred on 
the charge blobs or, more precisely, at the position of the 
density maximum of each blob. That is, we consider: 

V'^(r) = Aexp{-77i[(x - X,)^ + [y - Y,)^] - 772(2 - Z,)^} 

(10) 
where Ri = {Xi,Yi, Zi) gives the position of blob i, ini- 
tially set to coincide with the maximum of the corre- 
sponding electron density. The set of orbitals is orthogo- 
nalised using the Lowdin algorithm [33| , which preserves 
the equivalence of all orbitals, and then is normalised. 
The parameters 771 and 772 in the Gaussian exponent are 
varied in order to minimise the KS energy (at r^ = 40 
the optimal values of 771 and 772 are: 77ir^ = 2.3552, 
772r^ = 3.52). The approximation is remarkably accu- 
rate at rg =40, as confirmed by the low increase (less 
than 3 x 10~^ Ha per electron) of the optimal energy 
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(a) 
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(c) 

FIG. 7: (Color online) Magnetisation iso-surfaces for samples 
of iV = 32 electrons, (a) r^ = 25 (] m |= 0.3p6); (b) Vs = 35 
(I rn \— 2pb); (c) rs = 40 (| tti |= 2pb). Blue (dark) surface: 
spin up electrons. Yellow (light) surface: spin-down electrons. 
The Ts — 30 case is shown in Fig. [S] (c). 
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FIG. 8: (Color online) Density of states for the Kohn Sham- 
eigenvalues of fully spin polarised samples at different back- 
ground densities, A'^ = Nup = 32. The Fermi energy is ep, and 
wp is the plasmon frequency of the homogeneous electron gas 



at the same rs, given by ujp = \/3/' 



.3/2 



with respect to the plane wave estimate, based on the 
unconstrained optimisation of ^ 130000 coefficients per 
orbital. 

The Hessian d^Exs/dRfdRj for the iV = 24 supcrcell 
is computed by numerical differentiation of the energy 
upon moving the centre of Gaussians / and J by a small 
displacement S along the coordinates a and /3, respec- 
tively. The resulting matrix is combined with a kinetic 
part, to give the dynamical matrix of the chain, whose 
diagonalisation provides an estimate for the phonon fre- 
quencies. Analysis of the eigenvectors, or, more precisely, 
of the z-dependence of their polarisation vector, allows 
us to associate each eigenfrequency to a momentum q^, 
and thus to draw the dispersion relation over the first 
Brillouin zone. The results are shown in Fig. |9l Phonon 
frequencies belong to three branches, two of them being 
degenerate. The non-degenerate branch corresponds to 
vibrations along the z direction. It is an acoustic branch, 
whose frequency vanishes at g^ =0, then increases mono- 
tonically in moving towards the zone boundaries. The 
two other branches correspond to vibrations in the xy 
plane. The frequency of their q^ = modes does not 
vanish, because of the restoring force due to the external 
potential. Frequency, however, decreases with increasing 
\qz\, nearly vanishing at the zone boundary. This be- 
haviour clearly points to the easy deformation of the lin- 
ear chain obtained by displacing charge blobs perpendic- 
ularly to the wire axis, with nearest neighbouring atoms 
moving into opposite directions. The nearly soft mode 
thus corresponds to the formation of transversal dimers. 

This result suggests that the linear to non-linear tran- 
sition is second order, or at most weakly first order, in 




FIG. 9: (Color online) Phonon dispersion relation for a fully 
spin polarised sample of A^ = 24 electrons at Vs — 40. The 
ID Brillouin zone (along z) refers to the minimal unit cell of 
lattice constant a = L^/N, containing one single electron. 



qualitative agreement with the findings of Ref. [15| for 
the analogous transition in classical ion systems. 



IV. SUMMARY AND CONCLUSIONS 

Low dimensional assemblies of ions trapped into static 
and time-dependent electromagnetic fields have been ex- 
tensively investigated in the past, both computation- 
ally and experimentally, and could find applications in 
metrology and in quantum information. At the condi- 
tions of present experimerrts, the motion of ions is clas- 
sical, even though quantum mechanical aspects (strictly 
required for quantum information applications) are asso- 
ciated to the orientation and dynamics of spin moments. 

Recently, many experiments and calculations have 
been devoted, in particular, to analyse the so-called lin- 
ear to zig-zag transition [16| taking place in nearly ID ion 
systems upon varying the linear density of ions and/or 
changing the geometric parameters of the confining po- 
tential. In our paper, we have investigated the geo- 
metric and magnetic structure, electronic properties and 
low energy phonon-like excitations of nearly ID elec- 
tron systems confined by the electrostatic potential of 
a very thin cylindrical distribution of positive charge. 
We have considered a wide range of relative spin popula- 
tions iV„p and Ndn and focused on the low density regime 
(25 < rs < 40), where electrons tend to localise giving 
rise to blobs of negative charge distributed in space. In 
many respects, this model represents the quantunr coun- 
terpart of the classical ion systems. 

The results of our computations, carried within DFT- 
LSD, show that the quantum system displays a much 
wider variety of configurations and properties than in 



8 



the classical case. In fully spin-polarised systems, we ob- 
serve the stability of helicoidal density distributions for 
25 < r^ < 30, turning into a twisted string of localised 
charge dimers at lower density. Spin-compensated sam- 
ples at 25 < r^ < 30 display an intriguing double-helix 
structure, whose two strands have opposite spin polari- 
sation. The double helix unrolls into two linear chains of 
opposite spin at densities around r^ = 35, providing one 
of the few examples of virtually planar, zig-zag configura- 
tion found in our computations. Also in the case of glob- 
ally spin-compensated samples, charge dimers form at 
lower density (ts > 35) with a predominantly ferromag- 
netic coupling within each dimer and anti-ferromagnetic 
coupling among dimers. The spin-spin coupling, how- 
ever, is low at densities such that r^ > 30, and a siz- 
able amount of disorder is observed in the distribution 
of spins. At all densities, computations with a sin- 
gle spin-reversed impurity in an otherwise ferromagnetic 
sample reveal localised magnetic and structural defects, 
which might dominate the response of low density elec- 
tron chains to external perturbations. 

The density of states for the Kohn-Sham eigenvalues 
consists of several bands, some of them partially over- 
lapping, some other disjoint, according to density, aspect 
ratio and spin population. The ferromagnetic samples 
analysed in our study tend to be open-shell systems up 
to the lowest densities we investigated (r^ — 40), while 
spin-compensated samples are closed shell systems with 
a fairly wide gap separating occupied KS states from 
unoccupied ones. The open-shell character of the fer- 
romagnetic systems, together with the defects found in 
fully or partially spin-compensated cases, give rise to the 
non negligible low-frequency conductivity predicted by 
the Kubo-Greenwood formula for most of the samples 
investigated in our study. 

To the best of our knowledge, linear assemblies of 
electrons have never been made experimentally by the 
techniques used to trap atomic ions, although interest- 
ing experimental efforts on the single electron scenario 
are paving the way to an up-scaling. Our results pro- 
vide additional motivations for experimental investiga- 
tions along this line. Our findings suggest that electron 
systems of this kind, if ever realised, would exhibit a 



broad range of unusual and surprising properties, which 
could also find useful applications. 

From a fundamental standpoint, our analysis provides 
a useful complement to theoretical results obtained using 
classical simulation approaches by emphasising quantum 
mechanical and spin effects that might become impor- 
tant even in ion systems at sufficiently low T. Moreover, 
the density iso-surfaces computed in our study provide 
a more comprehensive view of the ground state proper- 
ties than the information given by low energy geometries 
of classical many-particle systems. The elongated struc- 
tures seen in some of the density iso-surfaces describe 
above, for instance, are a pictorial representation of low 
energy valleys in the potential energy surface, suggesting 
patterns for low frequency excitations. 

It is worth remarking that many-electron systems of 
this kind are already being made in condensed matter, 
using controlled doping to introduce mobile electrons into 
polymers and into semiconductor nanostructures j4l| . 
giving rise to quasi- ID structures whose size and den- 
sity are not far from those considered in our study. In 
these cases, however, the measurement of the properties 
of the system and their interpretation are certainly more 
difficult than in the ion-trap scenario. 

Our work paves the way to a variety of new stud- 
ies. For instance, an intriguing point to tackle would be 
the control of one of the "defects" found in our simula- 
tion samples. This could be realised by studying, via 
time-dependent density functional theory, the dynam- 
ics of the defect and its motion across the chain as in- 
duced by an external perturbative potential. It would be 
very interesting to study the effectiveness of this scenario 
for quantum communication protocols such as quantum 
state transfer, which are based on the use of quasi-lD 
lattices of interacting particles. 
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